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In this article we explore the phenomena of nonequilibrium stochastic process starting from the 
phenomenological Brownian motion. The essential points are described in terms of Einstein's theory 
of Brownian motion and then the theory extended to Langevin and Fokker-Planck formalism. Then 
the theory is applied to barrier crossing dynamics, popularly known as Kramers' theory of activated 
rate processes. The various regimes are discussed extensively and Smoluchowski equation is derived 
as a special case. Then we discuss some of the aspects of Master equation and two of its applications. 
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I. BROWNIAN MOTION : EINSTEIN'S THEORY 



A. Introduction 

The so-called Brownian motion was described for the first time in the year 1828 by the botanist Robert Brown. 
In investigating the pollen of different plants he observed that this become dispersed in water in a great number of 
small particles - the pollen grains. These were perceived to be in uninterrupted and irregular swarming motion. As 
the phenomenon repeated itself with all kinds of organic substances, he thought that he had found in these particles 
the 'primitive molecule' of living matter. He also found that all kinds of inorganic substances presented the same 
phenomenon and drew the conclusion that all matter was built up of 'primitive molecules'. 

Of the authors who carried out investigations on the Brownian movement before Einstein, we mention : Regnault 
(1858), Weiner (1863), Jevons (1870), Dancer (1870) and Delsaux (1877). 

The first precise investigation was due to Gouy (1888), who observed that the motion is more lively the smaller the 
viscosity of the liquid. He also ascribed the motion to the effect of thermal molecular motion of the liquid. Besides 
Gouy's work there was only other investigation of precise nature by Exner (1900) who showed that the velocity of the 
movement decreases with the size of the particle and increases with rise of temperature. 

Einstein was the first (1905) to formulate a correct picture of the entire problem. We discuss his theory in Sec. 1.2, 
1.3, 1.4. and 1.6. Q 

B. The irregular movement of particles suspended in a liquid and its relation to diffusion 

(a) Suppose there be suspended particles irregularly dispersed in a liquid. The particles are of microscopically 
visible size and their movements are of such magnitude that they can be observed under a microscope. We first 
consider here the irregular movements of particles which arise from thermal molecular movement. This gives rise to 
diffusion. 

Evidently it must be assumed (i) that each particle executes a movement which is independent of movement of all 
other particles, (ii) The movement of one and same particle after different intervals of time must be considered as 
mutually independent processes, if we think these intervals are not too small. 

We introduce a time interval r, which is very small compared to observed interval of time, but at the same time 
large such that the motion executed by a particle in two consecutive intervals r are mutually independent. 

Suppose that there are n suspended particles. In an interval r, x-coordinate of the single particle will increase by 
A. A has a different (positive or negative) value for each particle. For the value of A a certain probability law will 
hold. Let dn be the number of particles which experience a displacement between A and A + dA in time interval r. 
Then, 

nn 

— = 0(A) dA (1) 
n 

where the total probability is one, i.e., 

r+oo 



/+oo 
0(A) dA = 1 . (2) 
-oo 



Here 0(A) is the probability of jump of magnitude A for the particle, only differs from zero for very small values 
of A and fulfills the condition, 

0(A) = 0(-A) . (3) 

We confine ourselves to motion in one dimension (x). Let v = /(x, i), the number of particles per unit volume. We 
now calculate the distribution of particles at time t + r from a distribution at time t. 

We consider two planes perpendicular to x axis at x and x + A. Then the distribution at time t and space x + A 
evolves to a distribution at time t + r and at x as follows, 



1 We draw heavily from the two classic papers I and IV of Ref. jlj 
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/+oo 
f(x + A,t) 0(A) dA 
-oo 



(4) 



The integration over A takes into account of all possible jumps from x + A to x each with a probability 0(A). 
Since r is small we write 



f{x,t + r) = f{x,t) + 



df 
dt 



Again since A is small 

f(x + A,t) = f(x,t) + A 
Putting (5) and (6) in (4) we obtain 



df(x,t) A 2 8 2 f(x,t) 
dx 2! dx 2 



f+ m T 



f(x,t)cp(A) dA- 

p+oo 

f(x,t) 0(A) dA 



' f dj_ 

dx 

QfS f+oo 

dx 



A 0(A) dA- 
A 0(A) dA 



+ 00 



— OO 

:2 



dx 2 J 



dx 2 ) 2! 

+oo ^2 

■ rv-i • 



0(A) dA 
0(A) dA 



Since 



0(A) dA = 1 

-OO 



andQ 

i.e., A 0(A) = 0. We obtain 



0(A) = (—A) , 0(A) an even function , 



dt 



dj_ 

dt 



d 2 f 


r+°° a 2 


dx 2 „ 


/-oo ~V 




f+°° A 2 


}j 


"2l~ 



dx 2 



Putting 



We obtain the diffusion equation 



D = 



+oo 



A^ 
2! 



0(A) dA 



df(x,t) =D d 2 f(x,t) 



dt 



dx 2 



(5) 



(6) 



(7) 



(8) 



(9) 



where D is the diffusion coefficient. 

Note that both the equations (4) and (9) express the same law of evolution of distribution of particles in space and 
time. While (4) is of integral form, (9) is a differential equation. 

(b) Next problem that we investigate is how the distribution spreads in time. Mathematically speaking, this is an 
initial value problem. 

Let us suppose that at t — all the particles are concentrated at a point x = 0. This means that the density or 
the distribution function f(x,t) is infinite at this point and zero everywhere ( except this point ). We consider the 
problem of diffusion outward from this point. 



The jump of magnitude A has no preferential direction 
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At t = , f(x,0) =nS(x) (10) 

where by definition 

8(x) = , x ^ 
S(x) = oo , x = 

Also we have, 

/>+oo 
/(a;) 5 (a; - x Q ) dx . 
-OO 

Fourier representation of S(x) gives 
Writing f(x,t) as a Fourier transform, 



-00 



1 / , +°° 

/(*>*) =2^ / e ikx f(k,t)dk (11) 



and putting it in Eq.(9) we have 



2tt at 2tt 

Therefore, 

^M = -k 2 Df(k,t) (12) 

which can be solved to give 

f(k,t)=f(k,0)e- k2Dt . (13) 
Since f(x, 0) = n <5(x) and using the definition of 5 (or) we have 

n 5(x) = — / e lkx dk 
2tt J_ tx 

= — / e lkx n dk . (14) 

Again from (11) 



/(^,0) = — J e ikx f(k,0)dk 



hence by (10) we have 



Comparison with (14) gives 



1 



+OO 



nS(x) = / e lkx f(k,0) dk . (15) 

Z7T 



OO 



f(k,0) = n 



Therefore, 



•5 



f(k,t)=ne- k Dt . (16) 
Putting (16) in (11) we get the distribution at t [ see Appcndix-A for details ] 

The above equation shows how the distribution spreads in time (as shown in the figures ), 

(c) We now calculate the displacement X x in the direction of x-axis which a single particle experiences on an average 
- more accurately expressed - the square root of the arithmetic mean of the squares of the displacement, 



Thus (x 2 ) = 2Dt or 



1 f + °° 
(x 2 ) = - x 2 f(x, t) dx 

n J -co 

/+oo -1 
x 2 - 1 e- x ' 2 / iDt dx 

-oo V4tt Dt 

- / x 2 e- ax2 dx , a = l/ADt 

n J-co 



c + OO a 

—e- ax2 dx 



ad_ ■ " 

tt da 
a d 

w da V ol 

ir- ■ (18) 
2a v ' 



K = ^Ix 2 ) =V2~Wi . (19) 



The mean displacement is therefore proportional to the square root of time. This is a typical characteristic of 
Brownian motion. 



C. Diffusion and mobility 

We consider the suspended particles irregularly dispersed in a liquid. We consider this state of dynamic equilibrium, 
on the assumption that a force k acts on the particles and it depends only on position but not on time. For simplicity, 
we assume that the force is exerted everywhere in the direction of x-axis. 

We can look upon the dynamic equilibrium condition as a superposition of two processes acting in the opposite 
directions. 

(i) A movement of the suspended particle under the influence of the external force k. To be specific we consider 
a particle moving under the force of gravity k in a liquid. The particle is of spherical form with a radius p and the 
liquid has a coefficient of viscosity k. When the particle immersed in the liquid is falling under gravitational force k 
it experiences an opposing force of hydrodynamic origin. When the external force k balances the opposing force, the 
particle falls with a constant terminal velocity v . Thus if the external force of gravity is k and the opposing force is 
6ttkpvq, then 

^ = 7^~ ■ (20) 

blTKp 

If v is the no. of particles per unit volume, then v v number of particles pass a unit area per unit time ( or vk /(Sirup 
) under the action of external force k. 

(ii) The process of diffusion which is looked upon as a result of irregular movement of the particles produced by 
the thermal molecular movement of the liquid. 
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If D is the coefficient of diffusion of suspended and p the mass of a particle, then due to diffusion, the number of 
particles passing per unit area per unit time is 

Ox 

( — -D^r grams of particles crossing per unit area per unit time ) 



Under the condition of dynamic equilibrium we must have 

vk 



6lTKp 



ox 



At equilibrium, since the density under the force of gravity k varies as ( Boltzmann distribution ) 

k(x — xq)N 



v = Vq exp 



RT 



(21a) 



(21b) 



where N is the Avogadro number and R is the universal gas constant. Putting (21b) in (21a) we obtain, 

D=^^- (22) 
N Sirup y ' 

The coefficient of diffusion of the suspended particles therefore depends only on the coefficient of viscosity of the liquid 
and on the size of the suspended particles. 

If -^- p is denoted B ( "mobility" of the particle ) then (22) may be rewritten as, 



D=^B 

N 

This is the relation between mobility of the suspended particle and the diffusion coefficient. 



(23) 



D. Determination of Avogadro number 

We have found the diffusion coefficient D of a material suspended in a liquid in the form of small spheres of radius 
p as, 

D = ^ (^-) (24) 

Again the mean value of the displacement of the particle in x-direction in time t, 

\ x = V2Dt (25) 

By eliminating D we obtain 



r RT 1 , x 

This equation shows how depends on T, n and p. We will now calculate how great \ x is for one second. Take 
JV = 6x 10 23 

k = 1.35 x 10~ 2 poise (water at 17°c) ; 1 poise = lgm.cm _1 .sec _1 
p = 0.001 mm 

R = 8.31 x 10 7 erg. mole" Meg.KT 1 

We get 

X x = 8 x 10" 5 cm . 

On the other hand the relation (26) can be used for determination of N. We obtain thus 

X 2 X 3irnp 

where A^ has to be determined experimentally. 
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E. Experimental confirmation 



The first of the investigations confirming Einstein's formula in its original meaning was carried out by Seddig 
(1908) who look two photograph of an aqueous suspension of cinnabar on the same plate at an interval of 0.14 sec 
and measured the distance of the corresponding images on the plate. He found that on an average the displacement 
at different temperatures were inversely proportional to the viscosities as the theory demanded. 

Henri (1908) found similarly with the aid of cinematograph records of mean displacement of particles Caoutchonc 
that the time law, x 2 proportional to t, was followed. 

The establishment of first complete and absolute proof of the formula lies to the credit of Perrin and his group 
(1914), who followed the movements of single particle of gamboge or mastic under a microscope and recorded their 
positions at equal time intervals by means of an indicating apparatus. They determined N and found the value 
between 56 and 88 x 10 22 . 



F. Theoretical observations on Brownian motion and the existence of a random force 



It is now well established that irregular movements of the suspended particles in a liquid are caused by thermal 
motion of the molecules of the liquid. We now put forward two observations on the Brownian motion from a theoretical 
point of view to establish the existence of a random force. 

(i) From the molecular theory of heat we can calculate the mean value of the instantaneous velocity which the 
particle may have at the absolute temperature T. Thus the kinetic energy of the motion of a particle is independent 
of the size and nature of the particle and in depen dent of the nature of its environment, e.g., of the liquid in which the 
particle is suspended. The mean velocity \J (v 2 ) of the particle of mass m is therefore determined by the equation 

( y2 > 3 RT 

771——— = - — 28) 

2 2 N v ' 

with R = 8.3 x 10 7 erg.mole-^deg.K" 1 and N = 6 x 10 23 . 

We calculate (v 2 ) for particles in colloidal platinum solutions, for these particles we have the mass m = 2.5 x 10~ 15 
gm so that for T = 292 K 



J($*j= J— = 8.6 cm/sec (29) 
V mN 

(ii) We will now examine whether there is any prospect of actually observing this enormous velocity of a suspended 
particle. 

If we know nothing of the kinetic theory, we should expect the following thing to happen. 

Suppose that we impart to a particle suspended in a liquid certain velocity v by a force applied to it from outside. 
Then this velocity will die away rapidly on account of the friction of the liquid. The opposing force experienced by 
the particle is 6nnpv, where re =viscosity of the liquid, p =radius of the particle and v is the velocity of the particle. 
We obtain 

dv 

m— = —birnpv , (30) 



67rrep . 

v = v cxp[ 1\ . (31) 



On integration (30) gives 



m 

From this we calculate the time in which the velocity die away to one tenth of its original value, i.e., ii/io- From (31) 

v 6irnp 
In — = 1 1/10 

In 10 . , v 1 

h/w = ~ a 1— , with — = — 32) 

birKp/m v 10 

For platinum particles ( in water ) we have put 
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p = 2.5 x 1(T 6 cm 
k = 0.01 poise 
m = 2.5 x 10 -15 gm 



so that we get 



h/w = 3-3 x 10~ 7 



sec 



This means that the particle nearly completely looses its original velocity in the very short time ti/io through 
friction. But at the same time we must assume that the particle get new impulses from t he liq uid molecules during 
this time by some process that is the inverse of viscosity so that it retains a velocity y/ (v 2 ) on an average. But 
since we must assume that the direction and magnitude of these impulses are independent of the original velocity and 
direction of motion of the particles, we must conclude that the velocity and direction of the motion of the particle 
will be greatly altered in a very short time ti/ w and in a totally irregular manner. 

It is therefore impossible at least for ultramicroscopic particle to ascertain J (v 2 ) by observation. 

(hi) Although x/Jv 2 ) can not be observed, the change in position in time r ( which is much larger than ti/ w ) can 
be determined. We have already 



A* = vW^(V-) (33) 



N \37TK/9 

X x is the change in x-coordinatc in time r. Then the mean velocity in time interval r, we define as 



T 




(34) 



Since an observer can never perceive the actual path traversed in a arbitrary small time a certain mean velocity like 
^x/ T will appear to him as the instantaneous velocity. This is a measurable quantity. 

We have already seen that \J (v 2 ) is enormous for a particle when it is suspended in a liquid. On the other hand 
hydrodynamics tells us that this veloci ty if imparted to the particle by an impulsive force, will die away very rapidly. 
Therefore, to maintain the average \/ (v 2 ) the particle must experience impulses from the liquid molecules whose 
direction and magnitude are random. So to reconcile the kinetic theory with hydrodynamics we conclude that there 
must exist a random force F(t) acting on each particle, so that we write 

dv 

m— = -67TK/W + F(t) (35) 

Eq.(35) is an instantaneous description of motion of the particle. Here the average value of the random force is zero 

(F(t)) = and (F(t) 2 ) ^ 
Eq.(35) is the Langevin equation of motion for the particle. We discuss this equation in the next section. 

II. LANGEVIN DESCRIPTION OF BROWNIAN MOTION 
A. Introduction 

In the studies on Brownian motion we are principally concerned with the perpetual irregular motions exhibited 
by small grains or particles of colloidal size suspended in a liquid. As is now wellknown, we witness in Brownian 
movement the phenomenon of molecular agitation on a reduced scale by particles very large on a molecular scale - 
so large, in fact, as to be readily visible in an ultra-microscope. The perpetual motion of the Brownian particles is 
maintained by the fluctuations in the collisions with the molecules of the surrounding liquid. Under normal conditions 
in a liquid a particle will suffer as many as 10 21 collisions per second and this is so frequent that we cannot talk about 
separate collisions. And it is impossible to follow the path in any detail 

In the absence of any external force, one writes the Langevin equation for a free particle as 
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to— = —6nr)rv + F(t) 



(1) 



rj is the viscosity of the liquid, r is the radius of the particle and v the velocity of the particle. According to this 
equation the influence of the liquid medium on the motion of the particle can be split up into two parts : 

• A systematic part —6Trrjrv, represents a dynamical friction experienced by the Brownian particle. 

• A fluctuating force, F(t), which is characteristic of the Brownian motion. 
Two assumptions are made 

(i) F(t) is independent of v and (F(t)) = 0. 

(ii) F(t) varies extremely rapidly compared to variations of v. 

The second assumption implies that the time intervals of duration At exist such that during At the variations in 
v to be very small, while during the same interval F(t) may undergo many fluctuations. Although (F(t)) is zero, 
the average (F(t) 2 } does not vanish since negative swings of F(t) yield positive squared values. Suppose that the 
minimum time in which F(t) changes appreciably is called the correlation time, r c . The average of the product 
F(t)F(t'), vanishes for \t — t'\ > t c . Hence, (F(t)F(t')) , the correlation function of the random force is peaked about 
t = t' and falls off to zero in a time difference \t — 1'\ = t c . If r c is less than all other times of interest, e.g., 1/ (damping 
constant), we write 

(F(t)F(t')) = 2VS(t - t') (2) 

where V is some constant expressing the magnitude of the fluctuating forces. Eq.(2) along with (F(t)) = 0, completely 
defines the Langevin equation (1). 



B. General expression for mean square displacement 

The Langevin equation for the particle 

dv „ . 
to— = — mr\rv + F(t) 

gives — —6nrir(v} since (F(t)) = 0. 
Therefore, 



<«(*)> = M0)>exp(-^t) (3) 



TO 

Denoting Q-K-qr/m = T (damping constant) we rewrite Eq.(l) as 

mi = -mTx + F(t) . (4) 



since x = v, therefore, 



Multiplying both sides by x we get, 



x = -Tx+^- (5) 



TO 



F(t) 

xx — — Txx + x (6) 

TO 

Since x 2 = 2xx 
and x 2 = 2{x) 2 + 2xx 

Thus xx = i.x 2 — (x) 2 . 



With these relations we rewrite Eq.(6) as, 
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2 m 

or 



lx*-(i)* = -lj? + ±.xF(t) (7) 



x 2 - 2(±) 2 = -Pi 2 + —xF(t) (8) 



Taking the average we obtain, 



-^(z 2 > - 2((i) 2 ) = -r-^ 2 ) + -<zF(i)) (9) 
dr^ ar to 



Since[| 



i = w , 2((i) 2 ) = 2(v 2 ) = 2- 

m 

Also0 

= . 

Thus from (9) we have 

- ( , 2 > + r-(, 2 )-2^-.o. (to) 

We now solve Eq.(10) to obtain mean square displacement (a; 2 ). Put -^(x 2 ) = y and c = 2^^, then Eq.(10) gives, 

y + ry-c = (11) 



Let Ty — c = y' , then 



At t = 0, y' = A, hence 



From (12) 



y' + Ty = 

therefore y' — Aexp(— Tt) , A — constant . (12) 
Ty - c = A 

or A = -c ; c = 2(k B T/m) 



Ty-c = Aexp(-rt) 
c ^4 

y = p + p exp(-rt) 



Since, ^(a; 2 ) = y, 



On integration over to t we obtain, 



Since A = — c and c = 2^^, we obtain 



t+^-(l-e- rt ) (13) 



toF mr 2 

This is the general expression for mean square displacement of a Brownian particle suspended in a fluid 



4 



By equipartition of energy : |m(u 2 ) = ^fcgT, in one dimension 
The random force is independent of the position of the particle 
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1. long time limit 



When t — > oo, e rt = 0, we obtain the asymptotic behavior, 



2, 2fc B T j 



(x ) = — t (constant part is neglected) 



Putting T = we have, 



2x 2k B T : 



(x ) = — t • (14) 

The mean square displacement is proportional to time, t by (14). Comparison with Einstein's expression (x 2 ) = 2Dt 
gives 

&irqr 

which is the standard expression for diffusion coefficient calculated earlier by Einstein's method. 



2. short time limit 

When t — > small, from (13) 

2 2kBT 2kgT Tt 

we expand the exponential term to recover the leading order time dependence, 

(x 2 > = ^ 2 . (15) 

The mean square displacement in the short time is proportional to t 2 and is independent of the nature of the liquid. 
Or in other words, short dynamics is guided by inertial motion of the particle rather than any external influence. 



C. Relation between random and viscous force : The fluctuation-dissipation theorem 

We have the following Langevin equation 

6^™ + l_ 
at m m 

where, the first and the second terms are due to viscous and random forces, respectively. We rewrite (16) after 
multiplying both sides by v [ since T = (6Trrjr)/m ] 

v^- = -Tv 2 + —vF(t) (17) 
at m 



or 



Taking ensemble average we obtain 



^ - -2Tv 2 + -vF(t) 
at m 



l( v 2 ) = -2T(v 2 ) + ?-(vF(t)) . (18) 



The above equation requires the determination of (vF(t)); to this end we start from the identity 
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or 



t 

v(t') dt' = v(t) - v(t - At) (19) 

t-At 

v(t) = v(t - At) + [ v(t') dt' . (20) 

J t-At 

Multiplying both sides by F[t) and taking average 

(v(t)F{t)) = (v(t - At)F(t)) + f (v(t')F(t)) dt' . (21) 

Jt-At 

The first term of the right hand side of Eq.(21) vanishes because of the velocity v(t — At) at earlier instant t — At < t 
has no dependence on the future fluctuating force F(t). Thus, 

(v(t - At)F(t)) = . 

Therefore 

(v(t)F(t)) = f (v(t')F(t)) dt' . (22) 

Jt-At 



Putting the expression for v in above expression 



(v(t)F(t)) = 



t-At 
L 



1 



-Tv(t') + -F{t') 



m 



F(t) ) dt' 



T(v{t')F{t))dt'+ [ —(F(t')F(t)) dt' , (23) 

t-At Jt-At m 

the first term on the right hand side of Eq.(23) is zero again since t' is the earlier time (< t) and the fluctuating force 
at a later time t has no dependence on velocity earlier time t' except at t = t' for which the integral is zero. We are 
left with 

(v(t)F(t)) = - f (F(t')F(t)) dt' (24) 
m Jt-At 

We now assume that F(t), the fluctuating force, is stationary in time. This means the value of the integral (24) 
depends only on the difference t — t' but not on t and t' individually^ 



1 



(v(t)F(t)) = - / (F(t')F(t)) dt' 



t-At 
t+At 



= — / (F(t')F(t)) dt' 

2m Jt-At 

1 rt+At 

= Tr\ (F(i)F(t+s))ds 

2m Jt-At 

= ^ (F(0)F( S )) ds (25) 



t = t' + s 



Here it is important to note that, (i) since At ^> (t — t' = s), we put At ~ oo and (ii) The instant t is long ( equilibrium 
time ) and arbitrary so that we may put t = without any loss of generality. 

We now return to Eq.(18) and put the value of the (v(t)F(t)) from (25) to obtain 



5 In general, the stationarity of a correlation function means that it is invariant under time translation, i.e., (F(t)F(t')} — 
(F(t + T)F(t' + T)) 
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+oo 



!< W 2 )=-2r> 2 ) + l / <Is (F(U)F(s)) (2(i) 



In thermal equilibrium, the time derivative vanishes and the law of equipartition states 



l - m {J) = \k B T (27) 



Eq.(26) therefore gives, 



Applying (27) we obtain 



1 f +oc 

= -2I> 2 ) + — / ds (F(0)F(s)) 

m J — oo 



1 



+ 00 

ds (F(0)F(s)) . (28) 



2mfcsT J_ 

The above equation is known as the fluctuation-dissipation theorem, since it relates the dissipation T to the corre- 
lation of the fluctuating force F(t). It expresses a balance between the input of energy flow into the system ( particle 
) due to the fluctuating forces imparted by the liquid and the output of energy flow from the system due to dissipative 
forces exerted by the liquid on the system. 

One can easily have a relation between the diffusion in velocity space and viscosity. For this purpose we require 
that the fluctuating forces F(t) are instantaneously correlated, i.e., we use 

(F(0)F(s)) = 2V5(s) . 

Also since, T = (Qnr]r/m), we have from Eq.(28) 



6wr]r V 



/+oo 
S(s)ds 
-co 



m mk B T , 
or 

V = 6irr)rk B T . 

Thus this relation may be visualized as special form of fluctuation-dissipation theorem. 

III. BROWNIAN MOTION IN VELOCITY SPACE : FOKKER-PLANCK EQUATION 

A. Fokker-Planck equation 

In Einstein's method we considered the problem of Brownian motion in co-ordinate space, i.e., it concerns the time 
development of distribution of suspended particles in terms of f(x,t), or the probability of finding a particle at the 
position x at a time t. The law of evolution was stated to be, 

/+oo 
f(x + A,t)4>(A)dA (1) 
-oo 

where one takes into account of all the possible jumps of magnitude A from x + A to x, each with probability 0(A). 
The differential form of the above equation is the diffusion equation, 

df(x,t) d 2 f(x,t) 

~dT~~ u dx* (2) 

where D = J^°° ^-</>(A) dA is the diffusion coefficient in co-ordinate space. 

Herein we approach the problem of Brownian motion in velocity space (as in Langevin description ) and are 
concerned with the probability /(w, t) that a particle has a velocity v at a time t. The technique is applicable to any 
fluctuating quantity. For the sake of simplicity we consider, however, the problem in one dimension. 



14 



The time development of probability distribution f(v, t) of velocities may be stated as, 



f(v,t + r) 



f(v - A, t) 4>(v - A, A) dA 



(3) 



Here cf)(v — A, A) is the probability of a jump A for a particle with velocity v — A. Note that A has a dimension 
of velocity. Both r and A are small such that higher power of them may be neglected in the calculation. We now 
expand f(v, t + r) around t and f(v — A, t) (f>(v — A, A) around v such that 



(4) 



f(v - A, t)<Kv - A, A) = /(„, A) - A^P + A2 «W> 



9w 



<9v 2 



Putting (4) and (5) in (3) we obtain, 
/(»,i) + r|+0(r 2 ) 



J\v,t)<p(y, A) - A- 



9u 



2 <9v 2 

" + ■00 



dA 



/+oo q r /■+oo 

A) dA - — /(„, t) / A <f>(v, A) dA 

-oo L J — oo 

r /-+OC A 2 

/(«>*) y — <^(A,«)dA 



<9v 2 



(5) 



(6) 



Note that 4>(v, A) dA = 1 [ probability is normalized ]. We then write, 



1 f + °° (AM) 

- / A <p(v, A) dA = - — = Mi(«) 

T ./-oo T 



-OO 

1 /-+ 00 A 2 
- / — 4>{v,A)dA 



<A 2 ( V )) 



M 2 {v) 



Eq.(6) can be rewritten as, 



d.f(v,t) 
dt 



— M x (v) f(v,t) + —M 2 (v) f( Vl t) 



dv 



(7) 



M\(v) and M 2 (v) are called the drift and diffusion terms, respectively. Eq.(7) is called the Fokker-Planck equation. 
M 2 is the diffusion in velocity space and is not to be confused withD of Eq.(2) in co-ordinate space. 

For complete specification of Fokker-Planck equation one therefore needs the information about the drift, M\{v) 
and diffusion, M 2 (v) terms in Eq.(7). We calculate here these quantities for a specific model, i.e., the Brownian 
motion. 



B. Calculation of Mi(v) 



We integrate the Langevin equation 



- = -av + F(t) 



between t and t + t to obtain 



r-t+T 

v(t + r) - v(t) = -avr + J F(t')dt' 



(8) 



(9) 



a in Eq.(8) is the dissipation or damping constant and F(t) is the random force whose average is zero, i.e., (F(t)) = 0. 
We now put A = v(t + r) — v(t), therefore, 
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or 



is the drift term. 



(v(t + t) - v(t)) = (A) = -avr 



(A) 

Mx(v) = = -av 

T 



(10) 



From Eq.(9) we write 



C. Calculation of M 2 (v) 



A 2 = a 2 v 2 r 2 - 2avr J* +T F{t') dt' + G 2 t ( 



(11) 



where G t (r) = f t F{t') dt' . Since t is very small, the first term is neglected. Also the second term by virtue of 
(F(t)) = is zero. So we are left with ( after averaging ) 

(A 2 ) = (G 2 (r)) 



t + T r-t + T 



(A 2 ) = I dt 1 I (F(t')F(t")} dt" 



(12) 



If we assume that the random force is instantaneously correlated, i.e., 

(F(t')F(t")) = 2VS(t' - t" 

we obtain the average of A 2 as 



(A 2 ) = 2V / dt' / 5{t' - t") dt 



I J I \ JJ.II 



ft + T 

2V J dt' 



or 



(A 2 ) = 2Vt 



(13) 



(A 2 ) 

Therefore, M2(v) = = T>, the diffusion constant in velocity space, which is to be determined for the model 

represented by Eq.(8). 

To determine V we first integrate Eq.(8) formally with the integrating factor e~ at and obtain 



v(t) = v(0)e- at + e- at f e- at 'F(t') dt' 
Jo 

Rearranging and taking the square on both sides we have 

(v{t) - v{0)e- at ) 2 = e- 2at [ dt' f e a(t ' +t "^F(t')F(t")dt" 

Jo Jo 

Averaging yields 

((v{t) - v{0)e- at ) 2 ) =e~ 2at f dt' f e a(t ' +t "\F(t')F{t")) dt" 

Jo Jo 

Putting (F(t')F(t")) = 2V5(t' - t") in (16) we obtain 



(14) 



(15) 



(16) 
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((v(t) - v(0)e- at ) 2 } = 2e- 2at V f dt'e 2at ' 

Jo 

Explicit integration in the last equation yields 



(17) 



<(«(*) - v(0)e- at Y} = 2e~ 2at V 
V 



Jlat 



2a 



For 1/a, we have 



a 



V 



(1 - e- 2at ) 



(v(t) 2 } = - . 



a 



The above relation together with the equipartition theorem ^m(v 2 ) = \k B T leads us to 

V = ak B T/m 



(18) 



(19) 



Thus M 2 (v) = T> = akBT/m, ( diffusion coefficient in velocity space ). 

The above V is not to be confused D of Eq.(22) which is a diffusion coefficient for the Brownian particle in 
co-ordinate space. 

With these M\ and M 2 the Fokker-Planck equation for the Brownian motion can be rewritten as 



d.f(v,t) d d 2 



aknT 



m 



f(v,t) 



Putting p = mv the above equation (20) can be rewritten in the momentum space as follows 



dP(p, t) d 



dt 



d 2 

— [ap]P{p,t) + ^ 



dp 



dp 2 



aksT 



m 



P(p,t) 



(20) 



(21) 



where f(v,t) = P(p,t) represents the probability distribution function in momentum space. Here the underlying 
stochastic process is called the Ornstein-Uhlenbeck process (1930). 



IV. BROWNIAN MOTION IN PHASE SPACE (MOTION IN A FORCE FIELD) 

A. Kramers' equation 

We first derive here the equation of diffusion for an ensemble of particles with probability density distribution 
f(p, q, t) in phase space (i.e., q, p-space). The evolution of distribution from the time t to another time t + r is given 
by the following equation, 

/+oo 
f(p - A> q, t)(j){p — A, q, A)dA (1) 
-CO 

Had there been no Brownian motion, the motion would have been purely deterministic, i.e., 

I = $q) } (2) 

where JC(q) is the force acting on the particles. Thus the time development of q and p over a small time r would be 

qi = q + pr and Pi = p + JCt . (3) 

Here (q,p) and (qi,pi) specify the co-ordinate- momentum pair at time t and t + r. 

Because the particle is also subjected to a random Brownian force, we account for the all the possible jumps (of 
magnitude A) in momentum with the probability function <fr(p — A, q, A) and an integration over A in Eq.(l). 

Making use of Eq.(3) we rewrite Eq.(l) as, 
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/+oo 
/ (P - A, q, t)(f>(p - A, q, A)dA (4) 
-oo 



We now expand f(p + JCt, q + pr, t + r) in a Taylor series as, 



f{p , q ,t) + -)CT + -pT + -T + ... (5) 



Also expanding /(p — A, q, t)(f){p — A, q, A) as, 



/(P,9>*)W»9,A)- -^- A +2|-^2- A +•■■ ( 6 ) 



Integration over A in Eq.(6) gives 

e + OO 



f(p-A,q,t)4>(p-A,q,A)dA 

-OO 

/+oo o ^+oo o2 /-+00 a 2 

0(p, q, A)dA - — [/] y A<f>(p, q, A)dA + g^[f]j —HP, 9, A)dA 



/(p,5,i) - — \f{p, q,t)Mxb,q)\ + -^[f(p,q,t)M 2 (p,q)} . (7) 



We now put 



J^fog.AJdA = 1 

/+~A#p, g ,A)dA = A 

^^feq.A^A = A*/2 
Using Eq.(5) and Eq.(7) in Eq.(4) we obtain 

df df df 5 a ^ [J, ~ 

dp dq^ dt dp dp 

Dividing both sides of Eq.(8) by r we get 



_i/Cr + 4-pr + -^r = --[/A] + -j[/A*] (8) 



df df df d d 2 
where M x = A/t and M2 = A 2 /2r. Thus rearranging Eq.(9) we write 

Since the force /C(q) is derivable from a potential we write 

K(q) = -V'(q) (11) 
From the knowledge of Brownian motion (as calculated in the last section) we know, 



M% = — 7p = —77171; 
M2 = T> = mykBT 

Take the mass of the particle m = 1 for simplicity. Then writing q = x and p — v Eq.(lO) reduces to 



(12) 



The above equation is called the Kramers' equation ]6|,|7j. It describes the Brownian motion of a particle which is 
in a field of force. Or in other words a particle moves in an external field but in addition is subjected to irregular 
forces (at the same time) of the surrounding medium. While the first two terms are due to the deterministic motion, 
the third and the fourth terms are the drift and diffusion terms which are characteristic of Brownian motion. 
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B. Kramers equation as a generalization of Liouville equation and connection to equilibrium statistical 

mechanics 



The deterministic motion described by (13), i.e., 

corresponds to Liouville equation which forms the basis of equilibrium statistical mechanics under the condition 
-Jj: = ( which defines the equilibrium ) 

It is now wellknown that under equilibrium condition the distribution is a Maxwell-Boltzmann distribution, i.e., ( 
we assume m = 1 ) 

f(x,v) = Ze (15) 

where Z is the normalization constant. Thus (15) satisfies (14) for 5£ = as may be checked. 

It is important to emphasize that if we keep the Brownian dynamical terms in (13) as such and put the equilibrium 
condition ^£ = then it may also be checked that (15) satisfies (13) under the equilibrium condition ^ = . 

Kramers' equation may thus be regarded as a generalization of Liouville equation since it includes the Brownian 
motion in such a way that the basis of equilibrium statistical mechanics remains unaffected. 



C. Kramers' theory of activated processes 



Kramers' model for a chemical reaction consists of a classical particle of mass m (considered here to be unity) 
moving in a one-dimensional asymmetric double- well potential V(x). The particles co-ordinate x corresponds to the 
reaction co-ordinate and its value at the minima of the potential V(x), x a and x c denotes the reactant and the 
product states, respectively. The maximum of V(x) at x — xi, separating these states corresponds to the transition 
state (or activated complex) . All the remaining degrees of freedom of reactants and the solvent molecules constitute 
the surrounding medium whose total effect on the reacting particle is described by a fluctuating force and a linear 
damping. The correlation of fluctuating force gives rise to diffusion coefficient and the linear damping is responsible 
for the drift term. The stochastic dynamics for the reaction co-ordinate x and velocity v is governed by Kramers' 
equation n 



-ir v + ir{ v ( x +7^} 

ox ov 



d 2 



P(x,v,t) 



(1) 



The conditions are such that the particle is originally caught in the left well may escape in the course of time due 
to thermal activation by passing over the potential barrier. We want to calculate the probability of escape and its 
dependency on temperature and viscosity of the medium and compare the value with the result of 'transition state 
method'. The calculation rests on the equation of diffusion obeyed by a density-distribution of particles in phase 
space as written above [Eq. (1)]. 

To determine the steady state escape rate from A to C (say) we consider that there is a stationary situation in 
which a steady state probability current (flux) over B from A-^C is maintained. 

The stationary probability density must satisfy the following conditions : 
• Since we are considering a stationary situation, i.e., 4£ = we write, 



d_ 

dx 



— [V'{x)+ 1 v}+ 1 k B T-^ 



P{x,v) = 



(2) 



At the barrier top B we assume the linearized potential, i.e., we write by expanding V{x) around Xb 

(x - Xb) 2 



dV 

V{x) = V{x b ) + 



. 1 2 V 



here we have put f(x, v, t) = P(x, v, t) ; just a change in notation 
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Since 2£| = and - 

da; I a; = a; b Ox* 



io 2 we have, 



= V{x b ) - -u\{x - x h f 



(3) 



While considering the motion around x = Xb the above potential (3) has to be used in Eq.(2). 
• • • Near the bottom of the A-well all the particles are thermalized. Therefore we must have the usual Boltzmann 
distribution to be valid here, i.e. , 



P(x, v) — z 1 exp 



- l -V 2 + V{x)} n.;,.T 



0/t X X d 



(4) 



The linearization of potential has to be done at x = x a , i.e., we write 



dV 

V(x) = V(x a ) + ^~ 



. . i d 2 v 



(x - x a y 



Since §^ I = and — 

Ox I x=x„ 



= uj 2 a we have, 



V(x)=V(x a ) + -LJ 2 a (x-x a ) 2 . (5) 

While considering the motion around x — x a , the potential V(x) as given by Eq.(5) has to be used in Eq.(2) 
• • •• Near the bottom of the well C all the particles are (as if) removed. This implies the condition 

P(x,v)~0 for x > Xb ■ (6) 

Once the probability density P(x, v) fulfilling the above requirements is known, the population in the A-well n a 
and the flux j over the barrier will be given by 



//„ = ; dx dv P(x, v) 

I Awcll 
r+oo 

j = dv v P(x, v) 



Hence the steady state Kramers' rate is given by 

kA^c = j/n a . 
Our next task is to calculate j and n a separately. 



(7) 

(8) 

(9) 



1. Calculation of j 

Since we are considering the flux over the barrier B, the linearized potential V(x) — V(xb) — hui 2 (x — Xb) 2 has to 
be used. The Kramers' equation^] therefore reduces to 



d d 



v 7^ + ^7, i-^b( x ~ x b) + 7«} + 7^r— 



dx dv 

We now construct P{x, v) in the following form 

P(x, v) = v) exp 



P(x,v) = 



\v 2 + V(x) 



k B T 



X SS Xb 



(10) 



(11) 



7 since V'(x) = —u)\{x — Xb) 
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Putting Eq.(ll) in Eq.(lO) we get after a little bit of straightforward algebra 



d 



d 



if- 



Boundary condition for £(x, v) in Eq.(ll) should be such that 

• £(a;, v) — > 1 inside the well A, for 

• v) — > beyond the barrier top B, for x > Xb- 
We now use the following linear transformation 

u = v + a(x — Xb) 

where a is a constant to be determined later. This gives 



dv ' du 

Making use of Eq.(13) and Eq.(14) we obtain from Eq.(12) 



Z(x,v) =0 . 



or 



We now put : 



Since 



-av--[u J Ux-Xb) + ,v]- +1 k B T—=0 



r) 2 f Bf 



u)^(x — Xb) + v(a + 7) = —Am , A = constant (to be determined). 



u = v + a(x — Xb) 



we have, 



u>l{x - x b ) + v(a + 7) = —X[v + a(x - x b )} 



Comparing both sides 



From Eq.(17) we eliminate A to obtain 



which gives 



-A = a + 7 



a 2 + 7a — ojI = 



and 



Eq.(16) now reduces to 



A = —Lol/a 



d 2 £ d£ 



or 
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<9 2 £ Xu d£ 
— - H = 

du 2 jk B T du 

Our next task is to solve Eq.(20). We put |^ = y. Therefore Eq.(20) reduces to 

dy = 

du jk B T Uy 

Integrating over u Eq.(21) gives 

In j/ = — u 2 +lnF 2 , F 2 — constant of integration 

"fk B T 



y = F 2 exp 

Since 



Xu 2 



2 1 k B T 



( \u 2 

— =y = F 2 exp 



we get 

f u ( Xu 2 \ 
t;(u) = F 2 exp I — — - — — J + F\ , F\ = constant of integration 
Jo V 2 lk B l J 

We look for a solution that vanishes at x —* oo; the above integral should however remain finite for all |u| 
This implies A > 0, i.e., positive. Since 

A = —ujI/cl 

the negative root of a should be chosen to keep A positive, i.e., 



— W(i) 



2 



Thus A and a are known in terms of the given parameters 7 and co 2 of the problem. 
Next we determine F\ and F 2 (the integration constants). 

When x — > 00 then u — > —00, since w = u + a(x — x&) and a is negative. Again when x — > 00 we must 
— ► 0. Therefore from Eq.(22) we obtain 



/•- 00 / \ 2 
0-« + */ ^(-—1* 



or 

/•o 



Fl = ^4/_r exp (-^r) " u 

-Fi = F 2 > 



TT"fk B T 

2A 



Therefore we obtain 

= ^2 



ir-yk B T f u ( Xu 2 \ 



We then return to the expression for P(x, v) 
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P(x, v) = £(x, v) exp 
Using Eq.(26) we get from Eq.(27) 

P(x,v) =F 2 



\v 2 +V(x b )- \ul{x- x b ) 21 



2\ 



x exp 



1„,2 



k B T 



\u 2 \ 

6XP 



+ V(x b ) - \ujI{x - x h ) 
k B T 



, X ~ x b 



du 



Writing 



C 



njk B T 



and 



Eq.(28) reduces to 



F(x,v) 



2A 



rxl>( ~2^r) du 



P(xb,v) = F 2 exp 



V{x b ) 
k B T 



C exp --v 2 /k B T + F(x b ,v) exp --v 2 /k B T 



X « Xb 



The expression for the current j is given by 



J 



+00 



v P(xb, v) dv 



Since the first term in Eq.(29) can not contribute to Eq.(30) because of 

r +00 

ve-r v2/kBT dv = 



we write 



j = F 2 e k i 



+00 



F 2 e »bt (-k B T) 



v e~ v2/2kBT F(x b ,v) dv 



dv 



e -v /2k B T F ( Xb]W ) dv 



F 2 e~^ (~k B T) ^F(x,v)e- v ' 2/2kBT 
+°° OF 



dF 

dv 



-v*/2k B T dv 



F 2 e *bt (k B T) 



dv 



e -v*/2k B T dv 



Sincef] 



F(x,v) = l exp(-^) du 
dF _ ( \v 2 \ 

^ _cxp ^2^rJ ' 



(27) 



(28) 



(29) 



(30) 



(31) 



since at x — xt u — v 
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Therefore Eq.(28) gives 



j = F 2 e k B? k B T 



= F 2 e k B T 



+00 



/+00 
-00 



' 2-,k B T g 2k B T fa 



[2-,k B T + 2k B r] v 



dv 



= F 2 e~— {k B T) (2^k B T) 1 ' 2 ( 
Finally we get the expression for steady state current 

j = F 2 e-^B 1 (2^) 1 / 2 {k B Tf 2 ( 



1/2 



A + 7 



A + 7 



1/2 



Calculation of n a 



The number of particles in the left well A is given by 

^+00 p+cg 



/+00 r+00 
dv / dx P(x, v) 
-00 J — 00 



Since P(x, w) is given by 



P(x, w) = £(x, u) exp 
and £(x,v) is obtained from Eq.(26). Thus 



\v 2 + V{x) 



P{x,v) = F 2 
We have the following condition : 



/ irjk B T\ 



1/2 



k B T 



e 2~yk B T fa 



V 2A ; J0 

• As a; — > —00, i.e., the left well ; u — > 00 [since it = u + a(x — xt) and a is negative]. Therefore 



f 

Jo 



e 2 -y k BT du = 



2irk B T-y 



2A 



P(x,v) =F 2 



2it^k B T 



1/2 



k B T 



Since 



V{x) = V{x a ) + -u 2 a {x-x a y 



( ' 2Tt 1 k B T\ 1/2 _Vt£a) [ + Ca vL_ , /■ 

i a = f 2 f — j^—j e fcsT y e 2 " bt rfv y 



00 



F2 (^)" 2 e -^ (2 » 4i! T)./= (^) 



1/2 



= P 2 



(2irk B Tf/ 2 / 7 \i/2 



e k b J 
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We are now in a position to calculate the Kramers' rate 

k = j/n a 

From Eq.(34) and Eq.(36) we get 

F 2 {k B Tf 2 (2tt) 1 /2 

k = 



\l/2 v(*b) 



p (2*k B T)V> (1) l/2 ^ 



V(*a) 



(38) 



2?r \ a + 7 



1/2 



V(x b )-V(x a ) 
k R T 



(39) 



where F(x 6 ) — V(x a ) = E, the energy of activation. The pre-exponential factor in Eq.(35) can be simplified further 
as follows, 

Since 7 + A = a_ (negative root of a is o_) [see Eq.(17)] 

A = -( 7 + a_) 



7 + 



a + (positive root of a) 



A 



a+a+ 



A + 7 \a- J a_a + 
We thus get the final expression for Kramers' rate for arbitrary 7 



k 



2-KUJb 



7 
2 



(§)' 



-E/k B T 



We now consider the two limiting cases in the above equation : 
• when 7 — > (i.e., for small viscosity coefficient) 



k 



7^0 



<ff_ p -E/k B T 

2tt 



which is the transition state result (independent of 7). 
•• when 7 — > large, i.e., 7 » u)b (large viscosity limit) 



k 



7— > large 



27ro;5 



1 + IM 

+ 2 7 2 



-E/k B T 



LU a U) 



b e -E/k B T 



or 



27ro>b 7 



_ U^fc E /k B T 
n/-y_ > large — - """" ' 



2-7T7 



(40) 



(41) 



(42) 



The rate of reaction is inversely proportional to the viscosity. This observation has been corroborated by a number 
of experimental investigations in the recent past. 

The general result (39) which is valid in the intermediate to the strong damping limit provides a theoretical basis 
for Arrhenius expression for reaction rate k = Ae~ E / kBT proposed many years ago. 
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D. A simple connection to Transition State Theory 



We start from Kramers' equation which describes the motion of a particle in a force field governed by the potential 
V(x) simultaneously subjected to Brownian motion : 



d_ 

dt 



P(x,v,t) 



d d 



d 2 



-8- X V+ d-v {V/{x)+ ^ } + ^ BT dv 2 



P(x,v,t) 



We have shown that the equilibrium distribution P eq (x,v) corresponding to above equation is given by 



P eq (x, v) = Z e 



One assumes a that the particles initially residing in the left well are equilibrated and also that the above distribution 
is valid for all x ( around the bottom of the left well ) except at the barrier top x = B. We therefore put 

P eq (x,v) =0 for x > B 

The normalization constant Z is thus determined by the condition 

r +oo r B 



or, Zj_ 



+ 00 _1„2 

e k B 



/-(-OO P -D 

dv / P eq (x, v) dx = 1 
-00 J — 00 

~ dv j P eq (x,v) 
J —00 

, r B 

Z v/MbT / 

J — oc 



e k B T dx = 1 



e k B T dx = 1 



Expanding V^ir) around the left bottom x a 



V(x) = V(x a ) + -V"(x a )(x - x a ) 2 



we have [ write V"(x a ) = cj 2 ] 



Z y / 2irk B T e <=s T 



"bt dx = l 



or Z V^fcsT — e «bt =\ 



or Z 



[B -» +00] 
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Therefore the escape rate A; over the barrier is obtained by computing the outward flow over the top of the barrier. 



k 



/ v P(x b ,v) 
Jo 



dv 



Thus we have 



P{x b ,v) = Ze •'Bi- 



a.nd k as 



k = Z e 



k B T I y 

Jo 



e 2k B T dv 



2G 



Putting the value of Z in the above equation we get 

k = — e ~ Ea/kBT 

where E n — V(xb) — V(x a ) is the activation energy. 

This is the transition state result we derived earlier employing Kramers' method 7 = ( Note that this is not a 
dynamical theory like that of Kramers. So 7 does not appear in the theory. For 7^0, one has to consider the 
problem of energy diffusion. We state the main result; the rate constant k becomes proportional to 7 in this limit). 
This result implies that whenever the Brownian particle is at the top with positive velocity it will escape as if there is 
an absorbing wall at the barrier top. A rough interpretation of this transition state result is that the particle oscillates 
in an effective potential |w^(i — x a ) 2 provided by the left well and therefore hit the wall u a /2iT times per second and 
each time has a probability e~ E °/ kBT to cross over ^ 



V. OVERDAMPED MOTION : SMOLUCHOWSKI EQUATION AND DIFFUSION OVER A BARRIER 



A. Smoluchowski equation 



We wish to derive an equation of diffusion for an ensemble of particles with probability distribution function f(x,t), 
where the particles in addition to Brownian motion execute a deterministic motion in a force field. The potential is 
given by V(x). The evolution of distribution from time t to another time t + r is given by the following equation 



/+00 
f(x-A,t)<f>(A)dA . 



(1) 



In presence of Brownian motion, and under the influence of potential V(x) the equation of motion for the particle 
of unit mass (m = 1) is given by 



x + ~/x + V'(x) =F(t) 



where F(t) is the random force term. 

Under overdamped condition x <C jx we write 



jx 
or, x 



-V'(x)+F(t) 
V'{x) + F(t) 



(2) 

(3) 
(4) 



7 7 

The deterministic increment in x in time r corresponding to first term on the R.H.S. of Eq.(4) can be calculated as 



dx V'(x) 



(5) 



Expanding the functions / on both sides of Eq.(l) as usual we have 
,. , df df (dx 



+00 ■ 



r+00 r of 1 r) 2 f 

/+00 r+00 pp. r r+oo -i 

dA - J- J A 4>{A) dA + ^LJ_ -A 2 0(A) dA 



We now note 



1 



J^A^(A)dA = 
i J_+- A 2 0(A) dA = D) 



(6) 



(7) 



27 



Putting Eq.(5) in Eq.(6) and making use of Eq.(7) wo get 



d£ = _d£ ( V(x)\ | D d 2 f 

dt dx \ 7 / dx 2 

or 

df _ df fV'(x)\ d 2 f 

at {~r) +D dx^ (8) 

where D is the diffusion coefficient ( position ) and is given by k B T/-f as derived earlier ( Einstein's theory ). This 
equation is known as Smoluchowski equation. 



B. Diffusion of particles over the barrier 



We consider a particle moving in a potential field V(x) of the type shown in the Fig.(..). More generally, we may 
consider an ensemble of particles moving in the potential field V(x) without any mutual interference. We suppose 
that the particles are initially caught in the potential hole at x min . The general problem we wish to solve is the rate 
at which particles will escape over the potential barrier as a consequence of Brownian motion. 

The problem is very complex. However, considerable simplification can be made if we assume that the height of the 
potential barrier is large compared to the energy of the thermal motions, i.e., E ^> k B T. Under this circumstance, 
the problem can be treated in which the conditions are quasi- stationary. 

More specifically we may suppose that to a high degree of accuracy an equilibrium distribution exists in the 
neighborhood of x m i n . But this distribution is not valid for all values of x. We assume that beyond x max the density 
of particles is very small compared to the equilibrium value. And in consequence of this there will be slow diffusion of 
particles across x max tending to restore the equilibrium throughout. If the barrier were sufficiently high this diffusion 
will take place as though the stationarity prevailed. This condition is termed as a quasi-stationary condition. 

We thus consider the Smoluchowski equation, 



df(x,t) _ df 

dt dx 



V'(x)\ , k B T d 2 f(x,t) 



7 



7 



dx 2 



(9) 



Recasting the above equation in the form of a continuity equation we identify j as the current 

V'(x)f(x,t) k B T df(x,ty 



d£ = _d_ 
dt dx 
d 



7 



7 dx 



= -d-x 3{x > t) ■ 
In the stationary state j =constant, i.e., ^ = 0, where 



Rearranging the above equation as 



_ . _ V'{x)f(x,t) | k B T df 
7 7 dx 



df{x) + V'(x)f(x) _ 31 



dx k B T k B T 

and integrating between x min to A with the integrating factor 



(10) 



(11) 



(12) 



the equation (12) in the following form 



dx 



= e 



d_ 

dx 



k B T 



(13) 
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we obtain 



/(*) 



The constant current or flux 



31 
k B T 



f A +m 

/ e k B T 

J hin 



dx 



(14) 



/(or) e +fc s T 



f A e +7 ^ dx 



V(x) 



(15) 



Since /(x) at A is zero, i.e., f(A) = 0, we have 



J = 



I ^ mi ti ) 

k B T f(x mm ) e + *bt 



e 7 =b t dx 



(16) 



Around x m in, the current is almost zero. This defines an equilibrium condition in the neighborhood of x m in- Thus 
with j = the Smoluchowski equation yields ( see Eq.(ll) ) 



k B T df V'{x) 



or 



j dx 7 

df _ V'(x) 
dx k B T 

Integrating between x m i n to x ( a point in the left well ) 



In 



fix) 



fix) ■ 



V'jx) 
k B T 



dx 



or 



fix) fiXmin) e 



fx V'(x) 



dx 



-V(x) + V(x min ) 

fiXmin) e »B T 



The equilibrium population in the left well is given by 



f(x) dx = f(x mm ) 



where x\ and X2 are two points around x m i n . 
The rate of escape k is thus given by 



J Xl 



-V(x) + V(x min ) 



dx 



Thus from (16) and (19) we obtain 

k = 

















e k B T 




e >* b t 



k = j/n a . 



, y(% in ) 

fc R T 



- V(x) + V(x m j„) 

fc s T dx 



e +_ v" f e +fc s T dx f 2 e fc s T dx 



(17) 



(18) 



(19) 



(20) 



Therefore 
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We now make use of the following linearization of V(x) around x min and x max . For the integral 



we use 



and let 



dx 



V(x) = E - - wl (x - x max f 



-co and A — ► +00. Thus 

f +00 lEa-lcohx-Xnax) 2 ] 



r v 

J — 00 



dx = e 



e 2kBT '"t dx 



For the integral 



e k B T dx 



we use 



and let 



Therefore 



V ( x ) = o W ( X - x minf 



X\ — ► —00 
x 2 — ► +00 



e ! *u Tl dx = 



Putting the values of these integrals (23) and (25) in the expression for k in Eq.(21) we obtain 

Eg 

k B T e 



k 



7 



= ^0 U b e -E /k B T 
2 7T 7 

i) The rate of activation has a typical Arhenius form Ae~ Ea l kBT . 

ii) The rate is thus inversely proportional to the friction coefficient of the medium. 



(22) 



(23) 



(24) 



(25) 



(26) 



VI. THE MASTER EQUATION 
A. Master equation 

In statistical mechanics we deal with probability distribution functions. The master equation is a typical probability 
balance equation. 

Let us recall the good old theory of Brownian motion. Suppose a series of observations of the same Brownian 
particle gives a sequence of positions x\, x 2 , X3, .... 
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Each displacement Xk+i — Xk is an element of chance, i.e., is independent of earlier positions Xk—i, Xk-2, etc. This 
means probability distribution function does not depend on the previous history. Thus the position Xk+i depends 
only on xt- We call this stochastic process a Markov Process. Since many collision have already occurred during the 
displacement Xk+i — Xk, this displacement is much larger than the mean free path. 

We now recall the basic equation for evolution of probability ( Einstein ) distribution function f(x,t). 

f(x,t + r) = J f(x + A,t)cf>(A)dA (1) 

This equation relates the probability distribution function of a Brownian particle at x and time t + r, to that for 
the particle at a previous position x + A at an earlier time t. 4>(A) is the probability of a jump of magnitude A ( i.e., 
Xk+i — Xk ) in the picture from x + A to x. 

We now introduce the following notation for convenience. We denote 

y = x + A (2) 

and therefore, dy = dA. 
Also write 



0(A) = fa - y) , (3) 

where the arrow refers to the direction of jump from x — > y. Eqn.(l) can then be rewritten as [ denote f(x,t) by 
P(x,t)\ 

P(x, t + r) = j<f>{y^x) P(y, t) dy . (4) 
Expanding the left hand side as before we write 

P(x, t) + r^g^ =JHv^ x) P(y, t) dy . (5) 

Rearranging 



dt 

Now note that 

J 4>{x ^y)dy = J <X-A) dA 

= 1 (7) 



(Normalization of probability). 
Therefore we write 



dt 

Dividing both sides by r and rewriting 



we get 



dP(x,t) 
Jt 



(8) 



P(x,t)= J <t>{x^y) P(x,t) dy 
Putting (8) in (6) we get 

y ^ ' -'KJ 0( y _ X ) />,:,;. t) dy - j <t>{x - y) P(x, t) dy (9) 



W(y^x) = ^^ (10) 



= J W(y -» x) P(y, t)dy- J W(x -» y) P(x, t) dy (11) 
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W(y — > x) is the probability of a jump from y — > a; per unit time or the transition probability per unit time. 

The above equation is called the master equation. One can immediately write a discrete version of this equation as 
[ replace y by n, x by m and integral by summation ] 



dP m {t) 
dt 



In this form the meaning of master equation in very clear The first term is the gain of state m due to transitions 
from the other states n and the second term is the loss due to transitions from m to all other states n. Note that 
W nm > when n 7^ m. The master equation is thus a loss-gain equation for probabilities of separate states. 

The master equation is a doorway for studying the approach to equilibrium. The condition for equilibrium is defined 



Y^W nm P? = Y,W mn P% 



PP(i = m or n) must be identified with the equilibrium distribution function known from equilibrium statistical 
mechanics. 

The above condition states the fact that in equilibrium the sum of all transitions per unit time into any state m 
must be balanced by the sum of all transitions from m to all other states n. We now state another stronger condition 
that for each pair n, m separately the transitions must balance. 

W P eq = W P eq 

This is the principle of detailed balance and is true for all closed, isolated systems. 



B. Applications 

One step process : If we consider the jumps only between the nearest neighboring sites then the process is called 
an one-step process. 

The coefficient r n is the probability /time that being at n, a jump to n — 1 has occurred and g n is probability /time 
for a jump to n + 1 from n. For this the master equation (12) reduces to 

Pn = 9n-l Pn-1 + r n +l Pn+1 - (r n Pn + 9nPn) (13) 



1. Example 1 : unidirectional random walk 

Consider a one step process with constant transition probability. 

r n = , g n = q (14) 

The master equation is 

Pn = q{Pn-\ - Pn) (15) 

It is a random walk over integers n = 0, 1, 2, . . . with steps only to the right at random times. We start with a trial 
solution 



where a n (t) is to be determined. 
Hence putting (16) in (15) we obtain 



Pn = a n (t)e- qt . (16) 



a n = q a n -i 

= q a n - 2 

a>2 = q ai 

di = q a 
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Assume ao(t = 0) = 1 i.e., po(t = 0) = 1 as the initial condition imposed on (16). Therefore 

a\ = q t 

a 2 = q a\ = q 2 t 

which leads to 

a 2 = q 2 t 2 /2 
a n = q n t n /n\ 

Therefore the solution is 

n j-n 

Pn{t) = hr e ~ qt 

which is a Poisson distribution. Note that it is a non-stationary distribution. In the next example we consider a 
wcllknown stationary distribution. 

2. Example 2 : quantized harmonic oscillator interacting with a radiation field 

Consider n = 0, 1,2, . . . states of a harmonic oscillator. The oscillator is interacting with a radiation field. The 
interaction between the oscillator and the radiation field is causing the transition between the states of the oscillator. 
Energy of a state n is 

n+ 2~) hV ( ' 18 ^ 

Since the dipole moment matrix element between n — 1 and n is proportional to n, probability for a jump from 

n — 1 — > n 

9n-i =/3n (19) 

Probability of a jump from n — > n — 1 

r n = an (20) 
a, (3 are dependent on the frequency of light. We start from 

t>n = 9n-l Pn-l + r n +l Pn+1 - (ffn + T n ) p n 

p n = f3 n p n -i + a(n + 1) p n+1 - [(3(n + 1) + a n] p n (21) 
For a stationary distribution p n = 0. Therefore we have 

= f3 n p„_i + a(n + 1) p n+1 - [j3(n + 1) + a n] p n 

13 (n + 1) p n - a(n + 1) p n+1 = (3 n p n _i - a n p n = const = . (22) 

Hence 

(3 p n -! =ap n 

Therefore 

(3 

Pn = — Pn-l 

a 

which we rewrite as 
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p n = [ — ) Pa , Po = constant (23) 



a 



This distribution is an equilibrium distribution. We know from equilibrium statistical mechanics that 

Pn = const x e -"WfcT (24) 

Therefore equating (23) and (24) we obtain 



(e- hv ' kT ) 



or t = e -WW 



a 

Since g n is proportional to the radiation density p present 

= Cp 

For r n which is given by a there are spontaneous (A) and stimulated processes (Bp), 

a = Bp + A 

Therefore we have 



e -hv[kT 



A + Bp 
Cp 

Rearranging we write 

P = Ce -hu/kT _ B ( 25 ) 

which has the form of a Planck's distribution. A, B and C can be determined by comparing with Rayleigh-Jeans law 
and the Wein's law in long wavelength and short wavelength limits. 

APPENDIX A: EVALUATION OF EQ.(17) 

/+CO 
e- k2Dt dk 

/>+oo 

_ JZ_ p -x 2 /4Dt / e ikx e (kVDi) 2 e x 2 /iDt ^ 



2tt j.oo 

_ f,-x 2 /4Dt 
27T 



n 
2tt 



11 e -x 2 /4Dt 



- DC 
/■+OO 

2tt 



/ 

J — ( 

/ cxp [- UkVDt) 2 + (-ix/2VDt) 2 + (-ikx)\] dk 

/+oo 
exp -(k\/lTt-(ix/2\/Di)f dk 
-oo 

/+oo ^ 2 
-= e~ v dy ; put y = k\fWt - (ix/2\/Di), dy = y/Dt dk 
-oo y/Dt 



V ' V^Dt 

This result shows that a Fourier transform of a Gaussian function is a Gaussian. 
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